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Abstract 

We calculate the IP heavy quarkonium fine splittings at NLO and discuss the impact of the 
calculation on the Xf,(lP) splittings. 
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I. INTRODUCTION 



The calculation of the heavy quarkonium potentials and spectra in perturbation theory 
has a long history . Early calculations had to face two main problems: 

(1) the consistent inclusion of non-perturbative effects and (2) the poor convergence of the 
perturbative series. Recently remarkable progress has been made in both. Non-relativistic 
effective field theories of QCD, have provided a systematic way to factorize non-perturbative 
effects in heavy quarkonium observables. Moreover, they have proved useful in organizing 
and cancelling renormalon singularities and in resumming logs with renormalization group 
techniques: these being the main sources of large corrections in the perturbative series. For 
a recent review on effective field theories for heavy quarkonium we refer to Q]. 

This progress has triggered in the last few years a renewed interest in perturbative cal- 
culations of the heavy quarkonium spectrum. In (t| the bottomonium spectrum was 
calculated in a fully analytical NNLO calculation that implemented the leading order renor- 



malon cancellation. In 
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a numerical calculation of the spectrum was done that 



included NLO spin-dependent potentials. This calculation also relies (weakly, according to 
the authors) on some assumption about the long-range behaviour of the static potential. 

nn 

In |13l Il4j the hyperfine splittings of the IS" bottomonium, charmonium and be system 
were calculated at NLL accuracy. In general these calculations show rather stable (in the 
normalization scale) and convergent (in the perturbative series) results. 

In this work we calculate the n = 2, / = 1 quarkonium fine splittings at NLO. Some 
contributions were calculated in pL] , but a complete analytical calculation was missing up to 
now. The calculation may be relevant for heavy quarkonium states for which the momentum 
transfer scale p is much larger than the scale of non-perturbative physics Aqcd- Two situ- 
ations may occur under this condition Q|. Let us call E the typical kinetic energy of the 
heavy quark and antiquark in the centre-of-mass reference frame: in a non-relativistic bound 
state p ^> E. The first situation corresponds to quarkonium states for which E > Aqcd- 
Under this circumstance the heavy quarkonium potential is purely perturbative and non- 
perturbative contributions are parametrically suppressed with respect to the NLO correc- 
tions. The second situation corresponds to quarkonium states for which p ^> Aqcd E. 
Under this circumstance the potential contains a perturbative part and short-range non- 
perturbative contributions. A NLO calculation may help in this case to constrain the size of 
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the non-perturbative contributions affecting the spin- dependent potentials, which have been 
since long the object of intense study |16j]. Moreover, due to Poincare invariance, they are 
related by some exact relation to the non-perturbative contributions in the static potential 

121- 

The paper is organized as follows. In Sec.|n]we will calculate the fine splittings for n = 2, 
I = 1 quarkonium states at NLO. Our main result is Eq. (|26|). In Sec. IHII we will apply our 
result to the fine splittings of the Xf,(lP) states. 



II. n = 2, I = 1 FINE SPLITTINGS AT NLO 



We want to calculate the energy splitting 



(1) 



at NLO accuracy. In order to do this we need the relevant spin-dependent potentials and 
the wave function at NLO accuracy. 

The spin-orbit, Vls> an d tensor, Vt, potentials at NLO can be found, for instance, in jjjj]; 
they read: 



^Ls(r) 



V^(r)+V LS 
3C F a s (/i) 



2m 2 r 3 
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L S, 

y(lnr/i + 7 jB ) - ^(lnmr + 7 E ) - ^C A + |c F + ~2>n/ 



(2) 
(3) 

,(4) 



V T {v) 
Vi°\r) 

K«(r) 



^ 0) (r)+y«(r), 
C F a s (n) c 
4m 2 r 3 



'12, 



7T 



y(lnr/i + 7£) - C A (\sxmr + ~f E ) + ^C A + C F + ^T F n f 



(5) 
(6) 

(7) 



where L = r x p, S = Si + S 2 , S« = o-»/2, S 12 = 3r • a x r • <x 2 - cri • cr 2 , C A = N c = 3, 

C F = 4/3, Tp = 1/2, /?o = lliV c /3 — 4Tpn//3, n/ is the number of active massless flavours, 

a s is the strong coupling constant in the MS scheme and 7^ ~ 0.577216 is the Euler constant. 

In order to calculate the wave function at NLO we have to consider the heavy quarkonium 

r — I 

Hamiltonian at NLO. In the centre-of-mass reference frame it reads |l8j]: 



H 



m 



(8) 



C F a s (fi) 



6V®(r) = y (0) (r)^ 



7T 



y(lnr fi + lE ) + ^C A - ^T F n f 



(9) 
(10) 



We call \nljs) the eigenstates common to H, L 2 , J 2 ( J = L + S) and S 2 . At NLO \nljs) is 
given by 



\nljs) = \nljs} {0) + \nljs} {1 \ 



(11) 



where 



P!+y(°)(r)) \rdjs)® = E^\nljs)^, with = -m (%^Y , (12) 



For later use we notice that 



\nljs)W = J2 W s ') 

(kjtn),l'j's> 



,., t ^ Hkl'fs'\6V^(r)\rdj 8 )^ 



(13) 



(nljs\ f(r) L ■ S \n'l'j's') = S w S jf s J-^ + 1} /(/ ± 1} s(s + 1} (nljs\ f(r) \n'ljs)(U) 



(nljs\ f(r) S12 \n'lj's') = 5jfS ss '{Si 2 )ijs{nljs\ f(r) \n'ljs), 



(15) 



where 



(Si2>Zj S 



21+2 ■ 
21-1 ' J 



+2 



2/ 



21+3 ' 



/ - 1 

,3 = 1 
3=1 + 1 



, for I ^ 



(Sia)yo — (Si2/0is — ) 

and that 

<®(nljs\ f(r) \n'l'j's')V> = 5 lv 5 ]r 5j Q \nl\ f(r) \n'l)® , 

where \nl)(°> are the eigenstates common to [p 2 /m + V^(r)] and L 2 . 
The fine splitting E(l 3 Pj) - E(l 3 P f ) at NLO is given by 

£(1 3 P,) - E(l 3 P r ) = (21jl\V LS + V T \21jl) - (21fl\V LS + V T \21j'l), 



(16) 



(17) 



where 



(21il|V LSjX |21jl> = (°)<21il|VL S)T |21il>(°) + «<21il|Vig T |21jl>W + ^ <21il|^ T |21il> W. 

(18) 



Using the explicit expressions of the potentials, Eqs. (JU)-©, we obtain: 

(C F a s (/i)) 4 



(0) 



(21jl\V<$\21jl) 



v(0) 



m ■ 



256 



b'(j + i)-4] 



(°)(21il|y«|21il)(°) = (°)(21 J 1|^|21 J 1)(°) 

2/i 



7T 



*|ln. 

2 \ mCp(x s 



2 C 1 2 

3 A CpO-s 



3b S 



\T F n f 



(19) 
(20) 



(21jl|Kj 0) |21jl}(°) 



768 \ i5l2 / 1 J' 1 ' 



W(21jl|yW|21jl}(°) = (°)(21j'l|^ 0) |21jl)W 

2 



x 



7T 



2 V mCF<y s 



+ 1 -C A ln 



(21) 



(22) 



Equations (|2()jl and (|22j) correct the analogous expressions that may be found in 

To calculate ^ {21jl\V^ T \21jl)^ we use Eqs. (fTljl. (|T3jl . (fTB^l and the momentum-space 
representation given in 20]: 

(p|W) (0) (0) <W|q) (27r) 3 (5 3 (p-q) 



E, 



(0) 



(0) 



-7! /m — p 2 /m 
47rCV a s 



"7|/m — p 2 /m |p — q| 2 — 7f/m — q 2 /m 



1 They should read (for Z ^ 0): 

-m + 1) - Kl + 1) - .(. + 1)1 



x + 



2 V ttiCfol^ 



tp(n + ; + -0(2/ + 3) + V(2I) + 7£ 



n-Z-1/2 



V>(n + 2 + 1) + ip{2l + 3) + V(2Z) + 7£ - 



n-Z-1/2 



- — C A + -C F + -T F n f 
36 3 9 ' 



r/3 



x 1 



«s 
7T 



2 V mCpa s 



tp(n + Z + 1) + tp(2l + 3) + V(2I) + 1e 



1/2 



-CU In 



Cfol s 



ip(n + l + l)+ ip(2l + 3) + ^(21) + 7B 



1/2 



where ip is the derivative of the logarithm of the Gamma function. We thank F. J. Yndurain for commu- 
nications on this point. 
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1287rm 7 3 f 2 72 2 (p - q) 2 / 9 

(P 2 + 7l) 2 (q 2 + 7l) 2 I (P 2 + 7l)(q 2 + ll) V 2 



6 7 2 



2 6 7 2 



+ 



p 2 + 7 2 2 q 2 + 7 2 2 



3 4 7 2 4 7 2 
2 P 2 + 72 2 q 2 + 7 2 2 



with 



Co 



47 2 2 (p-q) 2 



-, and 7 2 = m C F a s /4. 



Eventually we obtain: 
'(21jl|^ s 0) |21jl) 



W/oi.-ilT/Wloi.-nW _ (o)/9i„-in/(?)|9i„-i\(i) - (°)(21jl|1/ L ( g ) |21jl) (0) 



x ■ 



'<21jl|V^|21jl> 

a 
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3 /155 2 2 . 2/x \ 31_ 

-A) 7T + In — f— H Ca 

4 P0 V 54 9 mC F aJ 24 



(p 2 + 7 2 2 )(q 2 + 7 2 2 ' 



\T F n f 



(23) 



(24) 



«(21jl|Kf 0) |21jl)(°) 



(°)(21jl|T4 0) |21jl) 



X ■ 



7T 



3^ /155 2 2 , 

~/3b vr + In 

4 M0 V 54 9 



CD = (°)(21j1|t4 0) |21j1>(°) 
2/x 



mC F a s y + 24^ " 6^ 

n 



(25) 



In Eqs. (I24J) and (125(1 the logs agree with what calculated in |4| using a variational method 
and the last two terms proportional to Ca and T F nj may be also derived from Eqs. 
(j!9|) and ()21|) . The other terms are the new contributions to the fine splittings at NLO 
provided by this work. 

Summing up all contributions, the fine splittings of the n = 2, I = 1 quarkonium states 
are given at NLO accuracy by 



P(1 3 P,) - E(l 3 P f 



m(C F a s 
768 



x<j3b'(j + l)-j"(/ + l)](l + ^ 



/3o 2 In. 



2/x 



+ 



173 TT 2 



mC F a s 36 3 

2 2 29 2 14 

— C A In + — Ca + -C F T F n f 

3 C F a s 18 3 9 ' 



+ [(S 12 ) lil -(S 12 ) lia ] 1 + 



7T 



A) 21b 



2/x 



mC F a s 



+ 



173 _ tt 2 
36 3 
4 



-Ca In ^ + ^Ca + C F - ~T F n, 



(26) 



Cj?a; s 6 

Equation (|2^|) is the main result presented here. 

The energy splitting between the centre-of-gravity energy P(l 3 P) c . . g . = [5P(1 3 P 2 ) + 
3P(l 3 Pi) + P(l 3 P )]/9 and the l l Pi state is entirely given at order (C F a s ) A a s /7i by the 



spin-spin potential at NLO and is, therefore, not affected by our analysis. It has been 
considered in jj, 0, 0] and we reproduce it here for completeness: 



£(l 3 P)c.o. g . - E(l 1 P 1 



m(Cp a s ) 4 a s 



288 



7T 



A) 7 ' 

T~1 Ca 



(27) 



III. APPLICATION TO THE x&(EP) SPLITTINGS 

In this section we apply Eq. (J26)) to the fine splittings of the bottomonium IP states, 
Xb(l-P). As argued in SecUJa necessary condition for this to make sense is that p ^> Aqcd- In 
order to have an idea of the size of the scale p, we may follow [4] and solve the (perturbative) 
self-consistency equation mCpa s {p) / (2n) = p for n = 2. We obtain p « 0.9 GeV for a choice 
of the mass of 4.73 GeV (i.e. half of the T(1S) mass). For comparison, in the charmonium 
ground state case we obtain p « 0.8 GeV for a choice of the mass of 1.55 GeV (i.e. half 
of the Jfij) mass). We see that it is reasonable to expect the typical momentum transfer 
scale for n = 2 bottomonium states to be larger than the scale of non-perturbative QCD 
( < 0.5 GeV). Noteworthy it is also larger than the typical momentum transfer scale of the 
charmonium^round state, which has been often assumed to be in the perturbative regime 
13j. Indeed, p has been assumed larger than Aqcd for Xi(lP) states in 




In Figs. □ and H we show P(l 3 Pi) - P(1 3 P ) and P(1 3 P 2 ) - P(l 3 Pi) respectively as a 
function of fi. The light and dark bands show the LO and NLO expectations respectively 
as from Eq. ()26|). The widths of the bands account for the uncertainty in a s (Mz) only. 
The mass has been chosen to be m = 4.73 GeV. The correct number of massless flavours 
in the bottomonium case is n/ = 3 We note that in the momentum region around 

the physically motivated momentum transfer of about 0.9 GeV the LO and NLO bands 
overlap with each other and with the experimental value. This means that the perturbative 
series shows convergence inside the physically motivated momentum transfer region and 
reproduces, inside the uncertainties, the experimental values. Note that the NLO curves 
have a relative maximum, i.e. a minimal sensitivity point, that goes from 0.82 GeV to 0.90 
GeVfor P(r 3 Pi)-P(l 3 P ) and from 0.83 GeV to 0.91 GeV for P(1 3 P 2 ) -P(l 3 Pi) as a s (M z ) 
goes from 0.1167 to 0.1207. Therefore, the physically motivated momentum transfer region 
overlaps also with the region covered by the minimal sensitivity points of the NLO curves. 

The agreement between the fine splittings calculated at NLO and the experimental data 
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FIG. 1: E(1 3 P\) — E(1 3 Pq) versus the normalization scale /x. The light band shows the LO ex- 
pectation, the dark one the NLO one as from Eq. \2b}) . The widths of the bands account for the 
uncertainty in a H (Mz) = 0.1187 ± 0.002 [21]. The horizontal band shows the experimental value 
32.8 ± 1.2 MeVlii]. 
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FIG. 2: E(1 3 P2) — E(1 3 P\) versus the normalization scale fi. The horizontal band shows the 
experimental value 19.9 ± 0.8 MeV j^l /. All the rest is like in Fig. Q 



appears reasonably good. Moreover it happens in a region where the physical momentum 
transfer is expected to be, the perturbative series shows convergence and the NLO curve is 
less dependent on the scale \i. However, the uncertainties are large so that no conclusive 
statement can be made at this point. A quantity less sensitive to ct s and to correlated 
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FIG. 3: p versus the normalization scale p. The horizontal line at 0.8 corresponds to the LO 
expectation, the curve to the NLO one. The band accounts for the uncertainty in aJMz) = 
0.1187 ± 0.002 The horizontal band at 0.607 ± 0.038 shows the experimental value l&l /- 

uncertainties in the splittings is the ratio 

_ E(1 3 P 2 ) - ff(l 3 P x ) 
9 ~ E{1 3 P 1 ) - P(1 3 P )' 8 

It has been suggested long ago that this observable may be rather sensitive to non- 

perturbative contributions in the spin-dependent potentials Q|. The ratio p is equal to 

0.8 at LO in perturbation theory, while it is measured to be 0.607 ± 0.038 for the IP bot- 

tomonium state j^]. The dark band in Fig. |3] shows the NLO expectation. We see that 

the NLO corrections go toward the data and may explain from about 15% to 65% of the 

difference between them and the LO value in the scale range 0.8 -j- 1 GeV. At the scale 

fj, = 0.9 GeV the NLO calculation seems to account for 20% -j- 40% of the difference between 

0.8 and the experimental value. At this stage we cannot say to what extent the remaining 

difference is due to non-perturbative corrections, since higher-order perturbative corrections 

not included in this analysis can still be large, as the scale dependence of p and the case of 

the n = 1 hyperfme splitting may indicate. In this respect it would be important to perform 

a NLL analysis similar to those done in [lsI . 14 1 . 
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